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Numerical and experimental data for the heat transfer as function of the Rossby number Ro in 
turbulent rotating Rayleigh-Benard convection are presented for Prandtl number Pr = 4.38 and 
Rayleigh number Ra = 2.91 x 10^ up to Ra — 4.52 x 10^. The aspect ratio V = D/L, where L is 
the height and D the diameter of the cylindrical sample, is varied between F = 0.5 and F = 2.0. 
Without rotation, where the aspect ratio influences the global large scale circulation, we see a small 
aspect-ratio dependence in the Nusselt number for Ra — 2.91 x 10^. However, for stronger rotation, 
i.e. 1/Ro ^ 1/Roc, the heat transport becomes independent of the aspect-ratio. We interpret 
this finding as follows: In the rotating regime the heat is mainly transported by vertically-aligned 
vortices. Since the vertically-aligned vortices are local, the aspect ratio has a negligible effect on the 
heat transport in the rotating regime. Indeed, a detailed analysis of vortex statistics shows that the 
fraction of the horizontal area that is covered by vortices is independent of the aspect ratio when 
1/Ro ^ 1/Roc- In agreement with the results of Weiss et al. (Phys. Rev. Lett., vol 105, 224501 
(2010)) we find a vortex- depleted area close to the sidewall. Here, we in addition show that there is 
also an area with enhanced vortex concentration next to the vortex-depleted edge region and that 
the absolute widths of both regions are independent of the aspect ratio. 

PACS numbers: 47.27.te, 47.32.Ef, 47.27.ek 
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I. INTRODUCTION 

The classical system to study turbulent thermal convection in confined space is the Rayleigh-Benard (RB) system, 
i.e., fluid between two parallel plates heated from below and cooled from above [Ull]- For given aspect ratio T = D/L 
{D is the sample diameter and L its height) and given geometry, its dynamics are determined by the Rayleigh number 
Ra = PgAL^/^tvu) and the Prandtl number Pr = v / Here ^ is the thermal expansion coefficient, g the gravitational 
acceleration, A the temperature difference between the plates, and v and n are the kinematic and thermal diffusivity, 
respectively. The case where the RB system is rotated around a vertical axis at an angular speed 1^, i.e. rotating 
Rayleigh-Benard (RRB) convection, is interesting for industrial applications and problems in geology, oceanography, 
climatology, and astronomy. The rotation rate of the system is non-dimensionalized in the form of the Rossby number 
Ro = j3g/S./ L/ {20)^ which represents the ratio between buoyancy and the Coriolis force. 

It is well know that the boundary layer behavior plays an important roll in the heat transfer properties of a RB 
system [Il[3]-[6]. However, several studies have shown that also the general flow structure in the system can influence 
the overall heat transfer. For non-rotating RB Sun, Xi and Xia [7 were the first to demonstrate that the global 
transport in any turbulent flow system can depend on internal flow structures. Later Xi and Xia [8] quantified the 
difference in Nusselt between the Single Roll State (SRS) and Double Roll State (DRS). Recently, Weiss and Ahlers [9] 
found that the difference in the heat transport between the SRS and DRS decreases with increasing Ra. Subsequently, 
Van der Poel et al. [10] found that in 2D RB convection a different flow structure can change the heat transfer up 
to 30%. An overview on the importance of the internal flow structure on the heat transfer is given by Xia [TP. The 
mechanism of the effect of the bulk flow on the heat transfer is as follows: More efficient bulk flow transports hot 
(cold) plumes to the top (bottom), thus creating a larger mean temperature gradient there, and thus a larger Nu. 
In this paper we will see that also in rotating RB convection the internal flow structures influence the overall heat 
transport in the system [121118]. 

It has been shown by several authors that three different regimes can be identified in RRB convection [T2"-T8^. As 
function of increasing rotation rate one first finds a regime without any heat transport enhancement at all in which 
the large scale circulation (LSC) is still present (regime I). Zhong & Ahlers [17] showed that - though the Nusselt 
number is unchanged in this regime - nonetheless various properties of the LSC do change with increasing rotation 
in this regime. Here we mention the increase in the temperature amplitude of the LSC, the LSC precession (also 
observed by Hart et al. [19] and Kunnen et al. [20 ), the decrease of the temperature gradient along the sidewall, and 
the increased frequency of cessations. The start of regime II (moderate rotation) is indicated by the onset of heat 
transport enhancement due to Ekman pumping as is discussed by Zhong et al. [13] and Weiss et al. [151 |2T]. When 
the rotation rate is increased in regime II the heat transfer increases further until one arrives at regime III (strong 
rotation), where the heat transfer starts to decrease. This decrease of the heat transfer in regime III is due to the 
suppression of the vertical velocity fluctuations [131 [III [22] . 

Several experimental [ij [HlTTl [23H26] and numerical [El Ulllll |20l|22l [ME] studies on RRB convection have 
shown that in regime II the heat transport with respect to the non-rotating case increases due to rotation. A detailed 
overview of the parameter ranges covered in the different experiments can be found in the Ra — Pr — Ro phase diagram 
shown in figure 1 of Stevens et al. [32 and the Ra — Ro — F phase diagram shown in figure [l] The heat transport 
enhancement in regime II is caused by Ekman pumping [13, 17^ 20, 23, 26] [271 [33]|34] . Namely, whenever a plume 
is formed, the converging radial fluid motion at the base of the plume (in the Ekman boundary layer (BL)) starts 
to swirl cyclonically, resulting in the formation of vertical vortex tubes. The rising plume induces stretching of the 
vertical vortex tube and hence additional vorticity is created. This leads to enhanced suction of hot fluid out of the 
local Ekman layer and thus increased heat transport. Corresponding phenomena occur at the upper boundary. 

Zhong et al. [T3^ and Stevens et al. [M] [22] used results from experiments and direct numerical simulations (DNS) 
in a F = 1 sample to study the influence of Ra and Pr on the effect of Ekman pumping. It was found that at fixed 
Ro the effect of Ekman pumping is largest, and thus the observed heat transport enhancement with respect to the 
non-rotating case highest, at an intermediate Prandtl number. At lower Pr the effect of Ekman pumping is reduced 
as more hot fluid that enters the vortices at the base spreads out in the middle of the sample due to the large thermal 
diffusivity of the fluid. At higher Pr the thermal BL becomes thinner with respect to the kinetic BL, where the base 
of the vortices is formed, and hence the temperature of the fluid that enters the vortices becomes lower. In addition, 
it was found that the effect of Ekman pumping is reduced for increasing Ra. This is because the turbulent viscosity 
increases with increasing i?a, which means that more heat spreads out in the middle of the sample. This explains 
why for the high i?a-number of Niemela et al. [35] no heat transfer enhancement was found in the rotating case. From 
this we can conclude that the internal flow structures influence the heat transfer in rotating RB as is also the case for 
non-rotating RB convection [THTOj. 

Recently, Weiss et al. [151 HI] showed that the rotation rate at which the onset of heat transport enhancement 
sets in (1/Roc) increases with decreasing aspect ratio due to finite size effect. This means that the aspect ratio is 
an important parameter in RRB convection. We report in this paper on a systematic study of the influence of the 
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□ DNS r=0.5 (Present paper) # EH r=1 .0-2.0 (Present paper) 

□ DNS r=1 (Present paper) # SB r=0.5 (Weiss and Ahlers (201 0)) 

□ DNS r=2 (Present paper) # SB r=1 .0 (Zhong et al. (201 & 2009)) 



FIG. 1. Phase diagram in Ra — Ro — F space for RRB convection with Pr = 4.38. The data points indicate where Nu has 
been experiment ahy measured or numerically calculated in a cylindrical sample with no-slip boundary conditions. The DNS 
data and the Eindhoven (EH) experiments are from this study. The experimental results from Santa Barbara (SB) are from 
Weiss et al. [15 and Zhong & Ahlers TT for the F = 0.5 and F — 1.0, respectively. The EH experimental data focus on F > 1, 
but for one case we also give the result for F = 1 for benchmarking with the earlier SB data, a) A three dimensional view on 
the phase space (see also the supplementary material), b) Projection on F — Ra phase space, c) Projection on 1/Ro — F phase 
space and, d) Projection on 1/Ro — Ra phase space. We also refer to figure 1 of Stevens et al. [32 for a full representation of 
the Ra — Pr — Ro phase space of RRB convection. 
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TABLE I. Simulations performed in this study. The columns from left to right indicate the following: the Rayleigh number 
Ra, the aspect ratio F, the inverse Rossby number 1/Ro, the number of Ro cases that is simulated (n^o), and the number of 
grid points in the azimuthal, radial, and axial directions (Ne x Nr x Nz) 
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aspect ratio on heat transfer (enhancement) for moderate and strong rotation rates, i.e., for 1/Ro ^ l/Roc We 
start this paper with a description of the experimental and numerical procedures that have been followed. Where the 
experimental and simulation results overlap an excellent agreement is found. Based on the numerical data we will show 
that for 1/Ro ^ XjRoc the Nusselt number is independent of the aspect ratio, while there are some visible differences 
for the non-rotating case. The reason for this is that in the non-rotating case there is a global flow organization, 
which can be influenced by the aspect ratio. In the rotating regime vertically aligned vortices, in which most of the 
heat transport takes place [36 , 37 , form the dominant feature of the flow. As this is a local effect the heat transport 
in this regime does not depend on the aspect ratio. In the last part of the paper we will analyze the vortex statistics 
for the different aspect ratios to support that the vortices are indeed a local phenomenon. 



II. EXPERIMENTAL PROCEDURE 

Various heat transport measurements on RRB convection were performed in the group of Guenter Ahlers in Santa 
Barbara. These measurements were done in an aspect ratio F = 1 sample and cover the Ra number range 3 x 10^ < 
Ra ^ 2 X 10^^, the Pr number range 3.0 < Pr < 6.4, and the XjRo number range < XjRo < 20. These data points 
are shown in the Ra — Ro — F phase diagram for RRB shown in figure [l] The experimental procedure that has been 
used in these experiments is described in detail by Zhong & Ahlers [17 . All the Nusselt number measurements of 
Zhong & Ahlers [17] were documented in the accompanying supplementary material of that paper, and we use those 
data in our figures here. Recently, Weiss & Ahlers [21 did similar measurements in an aspect ratio F = 0.5 sample. 
As these measurement data are not yet available in the open literature we do not include a comparison with their 
data. 

Zhong & Ahlers [TT restricted themselves to aspect ratio F = 1. Here we present new heat transport measurements 
from the Eindhoven RB setup [18 , which is based on the Santa Barbara design [l7l|3S], for aspect ratio F = 4/3 and 
F = 2.0. We also present one F = 1.0 measurement for one particular relatively small Ra number from the Eindhoven 
setup, (i) because this particular Ra number was not available in the Santa Barbara (SB) data set of Zhong & Ahlers 
[17] and we wanted to compare the results with those for F = 2 and F = 4/3 at the same Ra^ (ii) because we want to 
compare with the simulation results which are restricted to small Ra^ and (iii) because we wanted to benchmark the 
results of our new setup against those of the SB setup. 

The Eindhoven RB convection sample has a diameter D of 250 mm and due to the modular design of the setup 
measurements for different aspect ratios (F = 1.0 — 2.0) can be performed by replacing the sidewall. All measurements 
are performed with a mean temperature of 40°C (Pr = 4.38). For any given data point, measurements over typically 
the first four hours were discarded to avoid transients, and data taken over an additional period of at least another 
eight hours were averaged to get the average plate temperatures and the required power. Details about the setup and 
the experimental procedure, which is closely based on the Santa Barbara one [111 [38], can be found in Kunnen et al. 

m- 

Presently, the Eindhoven experiments can only be performed with one type of plate material (copper) and not with 
more as done in Santa Barbara. For this reason we can not apply the plate corrections to obtain the absolute Nusselt 
number without making some assumptions that cannot be verified at the moment. Therefore for the experimental 
results we restrict ourself to relative Nusselt numbers Nu{l/ Ro)/Nu{0). 

III. NUMERICAL PROCEDURE 

In the simulations the flow characteristics of RRB convection for Ra = 2.91 x 10^ — 4.52 x 10^, Pr = 4.38, 
< 1/Ro < 12.5, and 0.5 < F < 2.0, see also table [l] and figure [ij are obtained from solving the three-dimensional 
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Navier- Stokes equations within the Boussinesq approximation: 



Du 

DO 




1/2 



V^U + 6>£- -r—ZX U, 

no 



(1) 



{PrRay/' 



(2) 



with V • u = 0. Here £ is the unit vector pointing in the opposite direction to gravity, D/Dt = 9^ + u • V the 
material derivative, u the velocity vector (with no-slip boundary conditions at all walls), and the non-dimensional 
temperature, < ^ < 1. Finally, P is the reduced pressure (separated from its hydrostatic contribution, but containing 
the centripetal contributions): P = p — r^/(8i?o^), with r the distance to the rotation axis. The equations have been 
made non-dimensional by using, next to L and A, the free-fall velocity U = \^PgAL. A constant temperature 
boundary condition is used at the bottom and top plate and the side wall is adiabatic. Further details about the 
numerical procedure can be found in Refs. [39ti4T) . 

The resolutions used for the simulations are summarized in table |Tj These grids allow for a very good resolution of 
the small scales both inside the bulk of turbulence and in the BLs where the grid-point density has been enhanced. 
We checked this by calculating the convergence of the volume averaged kinetic €u and thermal cq dissipation rates 
as is proposed by Stevens et al. [42 . We find that these quantities always converge within a 5% margin at most. 
A comparison with the results of Stevens et al. [42] shows that this convergence rate is clearly sufficient to reliably 
calculate the heat transfer. As argued by Shishkina et al. 03] it is especially important to properly resolve the BLs. 
According to equation (42) of that paper the minimal number of nodes that should be placed in the thermal BL is 
Nbl = 4.7 {Nbl = 7.1) when Ra = 2.91 x 10^ {Ra = 4.52 x 10^) and Pr = 4.38. In the simulations at Ra = 2.91 x 10^ 
{Ra = 4.52 X 10^) we placed Nbl = 14 {Nbl = 22) points in the thermal BL, which is on the safe side. For the non- 
rotating case similar numbers are obtained for the grid point resolution in the kinetic BL. When rotation is applied 
the kinetic BL becomes thinner. This effect only becomes significant for the highest 1/Ro number cases considered 
here, and we emphasize that for all cases the number of points in the kinetic BL is above the criterion put forward by 
Shishkina et al. [43]. Furthermore, it is also very important to make sure that the results are statistically converged. 
Again we use the methods introduced by Stevens et al. [42 to check this. For all cases the convergence is in the order 
of 1% and it is much better for most. For the simulations at Ra = 2.91 x 10^ in the aspect ratios F = 0.5, F = 1.0, 
and F = 2.0 the average statistical uncertainty is only about 0.5%. 

As a last check we compare the numerical results with experimental data. For this we use the data of the group 
of Guenter Ahlers, who did high precision heat transport measurements in F = 0.5 [9|, F = 1.0 04j, and F = 2.0 [44] 
samples. For all cases we find that the numerical results are within 1% of these experimental data. For some cases 
the numerical data of the present simulations are even slightly below the experimental data. This is reassuring with 
respect to numerical resolution issues, as normally the heat transport in underresolved simulation is larger than the 
actual heat transport, see Stevens et al. [42] . 

In the simulations we partially neglect centrifugal forces, namely the density dependence of the centripetal forces, 
which in the Boussinesq equations show up as —2Fr rOr^ with the radial unit vector f [4^. Several authors [13! [2Ql[46] 
have shown that this is justified for small Fr numbers. For all experiments in the F = 1 sample, where the Fr number 
Fr = ^{Lj'i) I g is below 0.03 for all cases, this condition is fulfilled. For the experiments in the F = 2 sample 
the Froude number is below 0.05 for all experiments up to XjRo = 5. For higher 1/Ro the Froude number quickly 
increases up to 0.12 for Ra = 2.91 x 10^ and up to 0.49 for Ra = 1.16 x 10^. For the experiments at Ra = 2.91 x 10^, 
for which the highest Fr number is 0.12, this does not infiuence the results as we find a perfect agreement between 
the experimental and numerical results, see figure 2c. The experimental results for higher Ra in the F = 2 sample, 
where the Froude numbers are higher for the highest 1/Ro case, do not suggest a strong infiuence of this effect on the 
heat transport measurements. However, at the moment we can not rule it out completely either. 

We note that the simulations presented here are very CPU time intensive; about 1.5 million standard DEIS A CPU 
hours have been used, due to the large aspect ratios, the relatively high Ra numbers that are resolved with high 
resolution, the number of different Ro number cases, and the long averaging times that are needed to get sufficient 
statistical convergence. 



IV. RESULTS 



In figure [2] the heat transport enhancement with respect to the non-rotating case is shown for F = 0.5, F = 1, 
F = 4/3, and F = 2 and different Ra. Here we find an excellent agreement between experimental and numerical results. 
The figure shows that there is a strong heat transport enhancement due to Ekman pumping when 1/Ro > 1/Roc, 
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FIG. 2. The ratio Nu{l/Ro)/Nu{0) as function of 1/Ro for different Ra for a) T = 0.5, b) T = 1, c) T = 4/3, and d) 
r = 2. a) The results from the DNS at Ra = 2.91 x 10^, and Ra = 4.52 x 10^ are indicated by the black and magenta open 
squares, respectively, b) The experimental results for Ra — 2.99 x 10^, Ra — 5.63 x 10^ (run E4 of Zhong & Ahlers [17 ), and 
i^a = 1.13 X 10^ (run E5 of Zhong & Ahlers [17 ) are indicated in black, red, and blue solid circles, respectively. The DNS 
results for Ra — 2.91 x 10^, and Ra — 5.80 x 10^ are indicated by black and red open squares, respectively, c) The experimental 
results for Ra — 2.93 x 10^, Ra — 5.82 x 10^, and Ra — 1.16 x 10^ are indicated in black, red, and blue solid circles, respectively, 
d) The experimental results for Ra = 2.91 x 10^, Ra = 5.80 x 10^, and Ra = 1.16 x 10^ are indicated in black, red, and blue 
solid circles, respectively. The DNS results for Ra — 2.91 x 10^, and Ra — 5.80 x 10^ are indicated by black and red open 
squares, respectively. All presented data in this figure are for Pr — 4.38. 



where 1/Roc indicates the position of the onset of heat transport enhancement [HI [15]. The figure also shows that 
the heat transport enhancement decreases with increasing Ra. This is because the eddy thermal diffusivity is larger 
at higher Ra and this causes the warm (cold) fluid that enters the base of the vortices to spread out more quickly in 
the middle of the sample. This makes the effect of Ekman pumping smaller and this results in a lower heat transport 
enhancement. The breakdown of Nu at high 1/Ro is an effect of the suppression of vertical velocity fluctuations 
through the strong rotation. 

When the results for F = 0.5, F = 1, F = 4/3 and F = 2 are compared, they look pretty similar at first sight, i.e. 
a strong heat transport enhancement when 1/Ro > 1/Roc and a heat transport reduction for very strong rotation 
rates. However, when we look more closely we see that there are several important differences. First of all 1/Roc, the 
rotation rate needed to get heat transport enhancement, increases with decreasing aspect ratio. Weiss et al. [T5l [2T] 
already observed this and explained it by finite size effects. In addition, a close comparison between the F = 1 and 
F = 2 data in figure [2] shows that the heat transport enhancement is slightly larger in a F = 2 sample than in a F = 1 
sample. This can be observed in more detail in figure [3j where the relative heat transport enhancement at fixed Ra 
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FIG. 3. The ratio Nu{l/Ro)/Nu{0) as function of 1/Ro for Ra ^ 3 x 10^ and different T. The experimental results for T = 1, 
r = 4/3, and F = 2 are indicated in red, dark green, and blue solid circles, respectively. The DNS results for F = 0.5, F = 1, 
and F = 2 are indicated by black, red and blue open squares, respectively. All presented data are for Pr = 4.38. 
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FIG. 4. The absolute Nusselt number of the simulation data presented in figure |3] The non-rotating values are indicated by 
the filled symbols on the left-hand side. The results for F = 0.5, F = 1, and F = 2 are indicated by black, red and blue open 
squares, respectively. The black, red and blue dashed vertical lines indicate 1/Roc for F = 0.5, F = 1, and F = 2, respectively. 
All presented data are for Pr = 4.38. 



is compared for several aspect ratios. 

The absolute heat transfer shown in figure |4] is available only for the numerical case, (as explained in section 2). In 
that figure we see that there are some visible differences in the heat transport for the non-rotating case, i.e. without 
rotation we find that the heat transport in a F = 2.0 sample is approximately 5% lower than in a F = 1.0 sample. 
Similar differences in the heat transport as function of F have been shown by the numerical study of Bailon-Cuba et 
al. [47] and in experiments of Funfschilling et al. [44] and Sun et al. [48], although these experimental and numerical 
results seem to suggest that the heat transport becomes less dependent on the aspect ratio for higher Ra. The main 
point indicated by figure [i] is that the heat transport becomes independent of the aspect ratio once 1/Ro^ 1/Roc^ 
l/Roc ~ 0.14, 0.4 and 0.86 for F = 2, 1 and 0.5, respectively. We note that we also find for Ra = 5.80 x 10^ that the 
difference in Nusselt between the F = 1 and the F = 2 case is smaller for 1/Ro = 5 than for the non-rotating case, 
which is in agreement with the data presented in figure |4| 

We believe that the reason for this phenomenon lies in the flow structures that are formed. For the non-rotating 
case the flow organizes globally in the large scale convection roll. Because this global flow structure can depend on the 
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FIG. 5. 3D visualization of the temperature isosurfaces in the cyhndrical sample at 0.65 A (red) and 0.35 A (blue), respectively 
for Ra = 2.91 x 10^ Pr = 4.38, 1/Ro = 3.33 and F = 0.5 (left upper plot), T = 1.0 (right upper plot), and F = 2.0 (lower plot). 



aspect ratio, there can be small variations in the Nusselt number as function of the aspect ratio. For strong enough 
rotation, i.e. 1/Ro ^ 1/Roc, the global LSC is replaced by vertically- aligned vortices as the dominant feature of 
the flow [T3[ [m [161 nil Ell |46] . In this regime most of the heat transport takes place in vertically-aligned vortices 
j36l l37l l49l l5Q] . Because the vortices are a local effect the influence of the aspect ratio on the heat transport in the 
system should be negligible. This assumption is used in several models [^[371 [SI], which consider a horizontally 
periodic domain, that are developed to understand the heat transport in rotating turbulent convection. 

To investigate this idea we made three-dimensional visualizations of the temperature isosurfaces at Ra = 2.91 x 10^, 
Pr = 4.38, and 1/Ro = 3.33, for the different aspect ratios, see figure [s] Indeed the figures confirm that vertically- 
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FIG. 6. Visualization of the vertical velocity (left plot) and temperature (right plot) fields at the kinetic BL height for 
Ra — 2.91 X 10^, Pr = 4.38, XjRo — 5, and F = 2. Note that there is a strong correlation between the areas where the strongest 
vertical velocity and highest temperatures are found, namely in the vortices. Red and blue indicate upflowing (warm) and 
downflowing (cold) fluid in panel a (b). 




FIG. 7. The vortices at the edge of the kinetic BL as identified by the Q2D criterion for Ra — 2.91 x 10^, Pr = 4.38, XjRo — 5 
and F = 0.5 (upper left plot), F = 1.0 (lower left plot), and F = 2.0 (right plot). 
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FIG. 8. The dashed hnes indicate the average of the fraction of the horizontal area A covered by vortices for Ra — 2.91 x 10^, 
Pr — 4.38 and 5 < 1/Ro < 8.33 and different aspect ratios. Black, blue, and red indicate the data for F = 0.5, F = 1.0, and 
F = 2.0, respectively. The azimuthal average of the fraction of the area that is covered by vortices as a function of the radial 
position, i.e. {A)r/R, is indicated by the solid lines. Panel (a) shows the data as function of r/R and panel (b) as function 
of {R — r)/L close to the sidewall. For each aspect ratio the data are averaged for the simulation data at 1/Ro — 5 and 
1/Ro — 8.33. For each rotation rate 2 snapshots at the bottom and top BL are used. 



aligned vortices are formed in all aspect-ratio samples. Furthermore, the figure reveals that a larger number of vortices 
is formed in the F = 1 and F = 2 samples than in the F = 0.5 sample. This is expected, since these samples have a 
larger horizontal extension. 

In order to reveal the structure of the flow in more detail we show the temperature and vertical velocity fields at 
the kinetic BL height near the bottom plate for Ra = 2.91 x 10^, Pr = 4.38, 1/Ro = 5, and F = 2 in figure [o]. This 
figure clearly shows that hot fluid is captured in the up-going vortices. Similar plots (not shown) revealed that cold 
fluid near the top plate is captured in down going vortices. In order to study the vortex statistics one needs to have a 
clear criterion of what exactly constitutes a vortex. For this we use the so-called Q-criterion [151 iMl HSl IHl US] • This 
criterion requires that the quantity Q2D HHISl], which is a quadratic form of various velocity gradients, is calculated 
in a plane of fixed height. Here we always take the kinetic BL height. Following Weiss et al. [15] an area is identified 
as "vortex" when Q2D < —{\Q2d\)v^ where {\Q2d\)v is the volume- averaged value of the absolute values of Q2D' Here 
we have set the threshold for the vortex detection more restrictive than in [15 to make sure that only the strong 
upgoing (downgoing) vortices are detected. However, we note that similar results are obtained when a less restrictive 
threshold is used. 

The result of this procedure for Ra = 2.91 x 10^, 1/Ro = 5, and F = 0.5 — 2.0 is shown in figure [t] The figure shows 
that the vortices (both up and down going) are in general randomly distributed. However, note that no vortices are 
formed close to the sidewall. To quantify this we determined the radial distribution of the vortices from the plots 
shown in figure [7] and similar plots. The result is shown in figure [Sj Figure [8^ confirms that no vortices are formed 
close to the sidewall, while in the bulk their fraction is roughly constant. Figure [SfD shows that the size of the region 
close to the sidewall where no vortices are formed is roughly independent of the aspect ratio. This is in agreement 
with the predication of Weiss & Ahlers [21 that is derived from a phenomenological Ginzburg-Landau model. As 
detailed information about the flow field is needed to determine the vortex distribution it is very hard to obtain this 
data from experimental measurements [53 . In figure [sJd one can see that in the region {R — r)/L < 0.015 the value of 
(^4)^/^ decreases faster to zero than for {R — r)/L > 0.015, which is due to the vortex detection method employed in 
this study. More specifically, we detect only the core of the vortex. As the vortex core is always formed some distance 
away from the wall this causes an (artificial) enhanced decrease in the number of vortices that are detected in the 
direct vicinity of the wall. Due to the vortex detection method that is used it is also difficult to estimate the average 
radius of the vortices as only the core of the vortex is detected. 

However, and this is the main point here, one can see in figure [S] that the fraction of the horizontal area that is 
covered by vortices, see the dashed lines in the figure, is independent of the aspect ratio. This observation supports 
our finding that the heat transport is independent of the aspect ratio in the rotating regime. This result may seem 
somewhat unexpected based on the data shown in figure ^p. This figure namely shows that the absolute size of the 
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vortex depleted region close to the sidewall is approximately aspect ratio independent. Therefore one would have 
expected that the average horizontal area that is covered by vortices averaged over the whole area is higher for larger 
aspect ratio, because for larger aspect ratio samples this vortex-depleted sidewall region is relatively smaller than for 
smaller aspect ratio ones. However, just next to the vortex depletion region we find the vortex enhanced region. As is 
shown in figure [sJd, the absolute width of this region seems to be rather independent of the aspect ratio of the sample, 
too. Hence the effect of the vortex depletion and the vortex enhancement region on the horizontally averaged area 
that is covered by vortices cancel out in first order. 

In this paper we focused on the the influence of the aspect ratio on the fraction of the horizontal area that is 
covered with vortices. We note that Weiss et al. [15] have already investigated the effect of rotation rate on fraction 
of the horizontal area that is covered by vortices. There we showed that the average horizontal area that is covered by 
vortices increases approximately linear with 1/Ro when 1/Ro > l/Roc- This result is in agreement with the predictions 
obtained from a phenomenological Ginzburg-Landau like model that is discussed in that paper. Furthermore, it was 
shown by Kunnen et al. [53] that the vortex densities and mean vortex radii are mostly independent of the Taylor 
number Ta = Ra/{Ro^Pr) except very close to the bottom and top plates where more vortices are detected when the 
Taylor number is raised. As the numerical simulations considered here have been obtained for similar Ra and Pr as 
the ones used in the above studies we have not investigated the influence of these parameters on the vortex statistics. 



V. DISCUSSION 

In summary, we investigate the effect of the aspect ratio on the heat transport in turbulent rotating Rayleigh- 
Benard convection by results obtained from experiments and direct numerical simulations. We flnd that the heat 
transport in the rotating regime is independent of the aspect ratio, although there are some visible differences in the 
heat transport for the different aspect ratios in the non-rotating regime at Ra = 2.91 x 10^. This is because in the 
non-rotating regime the aspect ratio can influence the global flow structure. However, in the rotating regime most 
heat transport takes place in vertically- aligned vortices, which are a local effect. Based on the simulation results 
we flnd that the fraction of the horizontal area that is covered by the vortices is independent of the aspect ratio, 
conflrms that the vertically- aligned vortices are indeed a local effect. This supports the simulation results, which 
show that the heat transport becomes independent of the aspect ratio in the rotating regime. In addition, it conflrms 
the main assumption that is used in most models, which consider a horizontally periodic domain [MlETllM], that 
are developed to understand the heat transport in rotating turbulent convection. The analysis of the vortex statistics 
also revealed that the vortex concentration is reduced close to the sidewall, while the distribution is nearly uniform 
in the center. In between these two regions, there is a region of enhanced vortex concentration. The widths of both 
that region and the vortex-depleted region close to the sidewall are independent of the aspect ratio. This analysis 
highlights the value of numerical simulations in turbulence research: The determination of the vortex distribution 
requires detailed knowledge of the flow fleld and therefore it would have been very difficult to obtain this flnding 
purely from experimental measurements. 
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